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Abstract 

The electromagnetic form factors of the nucleon have been calculated 
in a chiral constituent-quark model. The nucleon wave functions are 
obtained by solving a Schrodinger-type equation for a semi-relativistic 
Hamiltonian with an effective interaction derived from the exchange 
of mesons belonging to the pseudoscalar octet and singlet and a lin- 
ear confinement potential. The charge-density current operator has 
been constructed consistently with the model Hamiltonian in order to 
preserve gauge invariance and to satisfy the continuity equation. 
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1 Introduction 

For a long time constituent- quark models (CQM) have been proposed to 
explain spectroscopic properties of hadrons within a nonrelativistic frame- 
work []]] (see also refs. 0, |3| and references therein). In these models the 
effective degrees of freedom are massive quarks moving in a long-range con- 
finement potential with the basic SU(6) spin-flavour symmetry. According 
to the analysis of ref. j|] the residual interaction responsible for the SU(6) 
symmetry breaking is described by the one-gluon exchange diagram and is 
identified with the hyperfine-like part of its nonrelativistic reduction. Given 
the fact that the mass of the three constituent quarks is small, the nonrela- 
tivistic approximation is not a priori justified. In addition some relativistic 
corrections are also necessary to account for the observed small spin-orbit 



effects |5], Relativized versions of CQM were then discussed [[7], |J. Al- 
ternatively, relativity is considered from the very beginning in the CQM 
adopting the light-front formalism H, ITU[ or using the socalled Bakamjian- 



Thomas construction to derive the Poincare invariant formulation of the 
quark model for baryons ||12|| . 

While rather successful in describing the octet and decuplet ground states, 
these models still face some problems, such as, e.g., the wrong level orderings 
of positive- and negative-parity excitations, which can be traced back to 
inadequate symmetry properties of the one-gluon-exchange interaction. 

The existence of an increasing number of near-parity doublets in the high- 
energy sector suggests that the approximate chiral symmetry of quantum 
chromodynamics (QCD) is realized in the hidden Nambu-Goldstone mode at 



low excitation and in the explicit Wigner-Weyl mode at high excitation [13 



Thus the spontaneous breaking of chiral symmetry and the associated ap- 
pearence of the octet of pseudoscalar mesons as the approximate Goldstone 
bosons induces a chiral interaction between quarks that is mediated by such 
mesons |14|| . Its spin and flavour dependence modifies the symmetry prop- 



erties of the Hamiltonian and ultimately leads to a correct ordering of the 



positive- and negative-parity states in the baryon spectra [15 . 

Various hybrid models have been constructed advocating meson exchanges 
in addition to sizeable contributions coming from gluon exchanges (see, e.g. 
refs. H,0,m,[Lp,gg). 

Recently, a chiral CQM has been proposed whose effective quark-quark 
interaction is derived from Goldstone-boson exchange alone involving the 
pseudoscalar meson octet and singlet |2l], [22]. The model is capable of pro- 



viding a unified description not only of the nucleon and Delta spectra but 
also of all strange baryons. 

A stringent test of the model would be to probe its eigensolutions in the 
description of the electromagnetic properties of baryons. In this paper the 
nucleon electromagnetic form factors are calculated without free parameters 
starting from the nucleon wave functions obtained with the model of refs. [fZI 



311 and using a charge- current density operator consistently derived along the 



lines proposed in ref. ||23|| . The model is briefly reviewed in section 2, while 
the expression of the charge- current density operator is given in sect. 3. The 
results are presented and discussed in sect. 4. 



2 The model 



The chiral model of refs. 21, 22] is semi-relativistic in the sense that the 



kinetic energy operator is taken in relativistic form: 



i=l 



with ?77.j the masses and Pi the three-momenta of the constituent quarks. 
This form ensures the average quark velocity to be lower than the light 
velocity, a requirement that is usually not fulfilled by nonrelativistic models. 
In addition, by the choice (|l|) one excludes negative-energy states ab initio 
and simply solves a Schrodinger equation for bound states without facing all 
the complications of a fully covariant treatment of the three-quark system. 
The dynamical part consists of a linear confinement potential, 

Vconiinj) = V + Cnj (2) 

depending on the interquark distance and the two fitting parameters Vo 
and C, and a sum of pseudoscalar meson-exchange potentials: 
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where <7j and Aj are the quark spin and flavour matrices, respectively. In the 
static approximation used in refs. [^T|, fffi] , the meson-exchange potentials are 
given by 

(5) 

with /x 7 being the meson masses and g 7 the meson-quark coupling constants 
(7 = 7r, K, rj, rj'). In the chiral limit there is only one coupling constant g$ for 
all Goldstone bosons. Due to the special character of the singlet rj' meson, 
its coupling constant go was allowed to deviate from g s . 
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Since one deals with structured particles (constituent quarks and mesons 
of finite extension, one has to smear out the 8 function in eq. ([|). In refs. |2T 
22fl a Yukawa-type smearing was used, i.e. 



g ^1 r ij 



47r S(fij) — > , (6) 



involving the cutoff parameters A 7 , which were assumed to follow a linear 
scaling with meson masses: 

A 7 = A + /t/i 7 . (7) 

Once the quark masses are fixed, the model has five fitting parameters, 
i.e. the depth Vq and the slope C of the confinement potential, the ratio 
do/ds of the singlet to octet meson-quark coupling constants and the two 
parameters Ao and k defining A 7 . 

The specific spin-flavour symmetry inherent in the chiral potential is re- 
sponsible for the correct level structure. According to the selected values of 
the quark masses different sets of numerical values for the five free parameters 
can produce fits to the baryon spectra of comparable quality. 

The Schrodinger-type equation for the model is accurately solved in the 
stochastic variational method P3 using wave functions that are expanded in 



basis functions involving correlated Gaussians as follows 

®jm,tm t (2k, Vk) = 4 vi exp{-(3x 2 k - 8yl + 7f fc • y k ) 
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depending on the Jacobi coordinates of partition k, 

Xk 

Vk = r k 



m p f p + m q r q 
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(9) 



where r*j and m 8 (i = 1, 2, 3) are the particle coordinates and masses, and 
{k,p,q) is an even permutation of (1,2,3). The represent the bipolar 
spherical harmonics 



LM L 



(10) 



The spin (isospin) parts arise from coupling single-particle spins (isospins) 
following the scheme 
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(pq) 

X(l/2,l/2)s mo 



X(l/2,l/2)* ^ 
Xl/2 <S> Xl/2 



'Xl/2 



SM S 



(11) 



For a given total angular momentum J and isospin T the stochastic varia- 
tional method selects basis functions according to a set of 6 discrete param- 
eters (L, A, I, S, so, to) and 3 continuous parameters 7, S). L is the total 
orbital angular momentum, A and I are the orbital angular momenta corre- 
sponding to Xk and y/~, S is the total spin, s the spin and to the isospin of 
the subsystem (pq) . 

The total wave function is composed of a symmetrized linear combination 
of basis wave functions of the form (||) . 

The Schrodinger equation can also be solved in momentum space. The 
basis functions are of the same form as in eq. (§) in terms of the corresponding 
Jacobi conjugate momenta {Px k ,Py k )- The method has the clear advantage 
of producing analytical solutions for the total baryon wave function both in 
momentum and space coordinates. 



3 The charge- current density operator 

The relativistic form of the kinetic energy does not permit the use of the 
traditional one-body current density operator, nor is it necessary to adopt so- 
phisticated procedures to include relativistic effects, such as those proposed, 
e.g., in refs. 0, ||, |25], |26], |27j. A gauge invariant charge- current density opera- 
tor can be derived consistently with the model Hamiltonian. It contains both 
one- and two-body terms. The one-body contribution includes the charge, 
the convective- and the spin-current operators. Their matrix elements be- 
tween free particle states are obtained following the functional derivative 
formalism proposed in ref. [p3|l - For a particle of charge e and mass m they 



are given in momentum space by the following expressions, respectively: 

m (x)\p) = e-^e^, (12) 



w&m = e ^(4 e '^' 2 ' (13) 

(j?, S '|J s (x)|p, S ) = _|_(y|* x (^-fl| s ) JLe^* (14) 



where 



Ep^Jp' + m?, E f = J$ 2 + m 2 . (15) 



With respect to the usual nonrelativistic expressions, in the semi-rela- 
tivistic approach only the spatial components of the charge- current density 
operator are affected, while the time component is simply given by the charge 
density. Furthermore, the energy denominator appearing in the current ma- 
trix elements reduces to 2m in the low-energy limit recovering the nonrela- 
tivistic approximation. 

The two-body current operator can be derived directly from the continuity 
equation consistently with the model Hamiltonian described in sect. 2 (see 



also ref. [[HJ). It is given by 



= -ieE{[- (1) xr( 2 )] 3 ^+(Ai 1) Ar ) -A«Af) S jK } 
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The A4 t 5 matrices are related to the SU(2) V-spin subgroup contained within 
SU(3). Together with the exchange character of the isospin dependence of 
the two-body current they result from the commutator between the quark- 
quark potential and the charge density that depends on A3 and Ag. Therefore 
only pion and kaon exchanges are allowed. 



4 The nucleon form factors 



As a first test of the model the electromagnetic form factors of the nucleon 
have been calculated. In this case only the one-pion exchange part of the 



two-body interaction, eq. (|3]), contributes. The calculation is fully consis- 
tent and without free parameters. In the following we give results for two 
parametrizations of the chiral constituent-quark model with different val- 
ues for the constituent-quark masses. The first parametrization is from 
refs. [^T[ with constituent-quark masses m u>d = 340 MeV. The second 



is a modified version with m u ^ = 250 MeV and accordingly readjusted pa- 
rameters of the model in order to obtain baryon spectra of similar quality. 

The electric (Ge) and magnetic (Gm) form factors are plotted in fig. 1 
for the proton and in fig. 2 for the neutron. The thin (thick) solid lines refer 
to m u4 = 250 (340) MeV. 

Two remarks can be made on the results. First, the fall-off of G P E and 
G^ 1 as a function of Q 2 is lower than observed. This Q 2 dependence reflects 
the fact that here the constituent quarks are assumed to be point-like. As in 
other constituent-quark models this assumption underestimates the electro- 
magnetic radii of the nucleon. Second, the value of G v ^ at Q = does not 
reproduce the nucleon magnetic moment. However, the ratio G^/G 1 ^ is in 
good agreement with the corresponding observed ratio of the proton to neu- 
tron magnetic moment, a feature common to all nonrelativistic constituent- 
quark models. The discrepancy at Q = is due to two effects, (a) Two-body 
currents do not contribute to the nucleon magnetic form factor because the 
Hamiltonian of eqs. ([3]) and (f|) does not contain the full axial dipole-dipole 
interaction that describes one-pion exchange completely. Therefore it is not 
possible to mix different values of the orbital angular momentum with the 
same parity, (b) The semi-relativistic form of the one-body current with 
an energy-dependent denominator suppresses its contribution with respect 
to the nonrelativistic case; it would therefore require a rather lower value 
for the constituent-quark mass in order to reproduce the nucleon magnetic 
moments. 

In order to improve the quality of the results without destroying the agree- 
ment with the observed baryon spectra one has to consider that constituent 
quarks are effective degrees of freedom with some spatial extension p8| , §]. 
As such, a charge form factor f(Q 2 ) could be appended to the charge density 
operator (|TJ) and the convective part of the current density operator (|TJ) as 
well as a magnetic form factor g(Q 2 ) to the spin part of the current density 
operator (0). This will modify the Q 2 dependence. In fact, a rather good 
agreement with data can already be obtained for G p ^ at Q 2 > 0.5 (GeV/c) 2 



assuming a simple dipole form factor 



= irn^p (17) 

common to ail(u and d) quarks. This is achieved with a rather small quark 
charge radius, i.e. r c = 0.35 fm, almost independently of its mass. Due to 
the small radius r c the Q 2 dependence of G P E , although largely improved, 
could not yet be reproduced. 

On the other hand, once constituent quarks are treated as extended ob- 
jects, it is not unreasonable to introduce an anomalous magnetic moment K 
in their electromagnetic form factor. Thus, besides a dipole form for f(Q 2 ), 
the following form for g(Q 2 ) has been considered: 

9(Q 2 ) = f(Q 2 ) + -j^Q^- (is) 

The actual value of k has been fixed in order to obtain the experimental value 
of the proton magnetic moment. For a quark mass m = 340 (250) MeV one 
obtains k = 0.867 (0.549). Correspondingly, the neutron magnetic moment 
turns out to be —1.828 (—1.812) n.m. in good agreement with experiment. 
The other two parameters a and b in eqs. (|17p and (|18D are then fixed by 
fitting the Q 2 dependence of G P M . The resulting values for the quark charge 
and magnetic radius are: r c = 0.691 fm and r m = 1.050 (0.935) fm with a 
quark mass m = 340 (250) MeV. It is worth noting that the extracted value 
of the quark charge radius is significantly close to the value required by the 
assumption of vector- meson dominance. Without any free parameter one can 
then calculate the other nucleon form factors. The results are shown in figs. 
1 and 2 by the thin (thick) dashed lines for m = 250 (340) MeV. A rather 
satisfactory agreement is obtained, especially for G E . 

The electric form factor of the neutron turns out to be too small in all 
cases. This is due to the deficiency of the charge density operator (|l|) de- 



rived in ref. [23]. In order to be consistent with the semi-relativistic Hamil- 
tonian no other contributions involving spin-dependent terms, like e.g. the 
Darwin-Foldy correction, are possible and the quark charges add up to a 
total vanishing neutron charge. 



5 Conclusions 



A completely consistent calculation of the nucleon electromagnetic form fac- 
tors has been performed within the chiral constituent-quark model proposed 



in refs. [pi] , |22| . Considering point-like quarks is not sufficient to reproduce 
the observed form factors. In particular, the model eigenfunctions do not 
permit contributions from two-body currents arising from pion exchanges 
and the semi-relativistic one-body current alone fails to produce the correct 
values of the nucleon magnetic moments. However, with the inclusion of 
suitable electromagnetic form factors for quarks and considering an anoma- 
lous quark magnetic moment a rather satisfactory agreement with data is 
obtained. 

Possible improvements of the dynamic model are obviously under discus- 
sion. The pseudoscalar tensor term has been neglected in the model Hamil- 
tonian of refs. f2l| , |22j . In ref. |29] its effect on the light-baryon spectrum 
was estimated together with that of the Thomas-Fermi precession spin-orbit 
contribution arising from the scalar confining interaction. The result was 
that the agreement with the observed spectra was destroyed. The inclu- 
sion of vector- (p, u>, <fi, K*) and scalar-meson (a) exchanges was considered 
in ref. [BO]. The tensor forces of vector- and pseudoscalar-meson-exchange 
interactions have opposite signs and largely cancel each other. The effects of 
the spin-orbit forces from the vector- and scalar-meson-exchange interactions 
are rather weak. The problem of the spin-orbit force from the Thomas- Fermi 
precession (which was not taken into account in ref. pOfl ) remains, however. 
In any case, pseudoscalar tensor and vector-meson exchange contributions 
modify the model eigenfunctions, so that one can expect that the two-body 
currents will contribute even in the simplest case of the nucleon electromag- 
netic form factors. Work along these lines is in progress. 
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Figure 1: The electric (G P E ) and magnetic (G P M ) form factors of the proton 
as a function of the four- momentum squared Q 2 . The thin (thick) solid lines 
refer to a quark mass m u ^ = 250 (340) MeV. The thin (thick) dashed lines 
for m u ^ = 250 (340) MeV include the effects of electromagnetic form factors 
for quarks (see text). Experimental points are from ref. j31| (solid circles), 
ref. |32| (open circles) and ref. |$3[ (triangles). 
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Figure 2: The same as in fig. 1, but for the neutron. Experimental points for 
G% are from ref. |35| (open circles), ref. |36[ (solid circles), ref. |37| (cross), 
ref. |38| (triangle) and for G\,j from ref. (crosses), ref. |33| (triangles), 
ref. p5| (open circles), ref. |3(J (solid circles), respectively. 



